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Abstract 

We use functional renormalization group methods to study gravity minimally coupled to a free 
scalar field. This setup provides the prototype of a gravitational theory which is perturbatively non- 
renormalizable at one-loop level, but may possess a non-trivial renormalization group fixed point 
controlling its UV behavior. We show that such a fixed point indeed exists within the truncations 
considered, lending strong support to the conjectured asymptotic safety of the theory. In particular, 
we demonstrate that the counterterms responsible for its perturbative non-renormalizability have 
no qualitative effect on this feature. 
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I. INTRODUCTION 



Quantized general relativity is notoriously non-renormalizable at the perturbative level. 
Such an understanding has been achieved after a number of celebrated calculations that, 
starting with 't Hooft and Veltman's seminal work |l|], have disclosed the appearance of 
non-renormalizable divergences already at one- loop in the presence of matter l|, |2j , and at 
two- loop for pure gravity 3|, |4|. The situation is neither improved by the presence of a 
cosmological constant [5|, nor by non-minimal couplings 

aa. 

The general conclusion usually taken out of these results is that general relativity is not 
'undamental and can only be quantized as an effective field theory. In this approach (see 
8, the gravitational action is organized in an energy expansion in curvature invariants. 
Once the scale for an experiment is identified, only the pertinent terms are then retained, 
allowing one to make predictions. A problem arises, however, once the energy is such that 
the curvature in Planck units reaches unity. At this point all curvature invariants are of the 
same order and an infinite number of couplings has to be fixed, so that the predictive power 
is lost. 

A different conclusion can be attained if instead gravity turns out to be asymptoti- 
cally safe (AS) [lo| (see [ll, 12, [k], [ijj] for recent reviews). This scenario is based on 



Wilson's modern viewpoint on renormalization [15| and envisages the existence of a non- 
Gaussian fixed point (NGFP) of the renormalization group (RG) flow with a finite number 
of ultraviolet-attractive (relevant) directions. For RG trajectories attracted to the NGFP 
in the UV (spanning the UV critical surface of the fixed point), the fixed point ensures 
that the theory is free from uncontrollable UV-divergences, while the finite dimensionality 
of the surface ensures the predictivity of the theory at all energy scales. These criteria rep- 
resent a non-perturbative analogue of the requirements underlying the usual perturbative 
renormalizability, which is recovered in the case of the fixed point being the Gaussian one. 
In recent years, significant evidence for the asymp totic safety of gravity has been gath ered 
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28, 



29|, 



by use of functional RG techniques 
though support also comes from lattice simulations 30j. The former approach generally 
employs a Functional Renormalization Group Equation (FRGE for the effective average 
action originally derived in [3l| and first applied to gravity in 16j. Since an analysis based 
on the full equation is probably impossible, investigations usually rely on truncations of the 
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theory space, whereby only a finite number of interaction-terms are retained. The reliability 
of the results found within such approximations can then be supported by considering their 
stability under a gradual extension of the truncation subspace. Indeed, all truncations 
studied so far, from the Einstein-Hilbert-, to the R 2 - and general f{R)-, up to the R 2 + C 2 - 
truncations, give rise to a coherent picture, pointing at the existence of a NGFP dominating 
the UV behavior of gravity. 

A possible criticism on these results is that they are based on truncations which only 
contain interactions that are also unproblematic for the on-shell perturbative renormaliz- 
ability. It is therefore a fundamental test for AS to include potentially dangerous terms in 
the truncation ansatz and study their effect on the fixed point structure of the theory. In 
pure gravity, the first non-trivial counterterm would be the Riemann-cube term of |3j, ]4\. 
Including this term in the truncation ansatz is, however, technically very involved and be- 
yond the current FRGE-techniques, even though the work presented in [28j|, which for the 
first time permitted us to go beyond the class of f(R) -truncations, constitutes significant 
progress in that direction. 

A technically less demanding, but equally illuminating, alternative is to study trunca- 
tions for matter-coupled gravity. In this case the non-renormalizable counterterms already 
appear at one-loop and the occurrence of divergences proportional to R 2 and C iiVpa C ilvpa , 
which do not vanish on-shell, signal the break down of perturbative renormalizability. To 
date, investigations of matter-coupled truncations, while also corroborating the asymptotic 



safety scenario, have remained restricted to the Einstein-Hilbert case [27I. l32l. l33l. I34J. l35| . In 
the present paper, we go beyond this restriction, and study the non-perturbative RG flow of 
gravitational higher-derivative terms in the presence of a free, massless, minimally coupled 
scalar field, cf. eqs. ffl2l and (1251) below. Anticipating our main result, we find that the 
NGFP previously reported for the Einstein-Hilbert case persists under the extension of the 
truncation subspace. This constitutes further evidence for the non-perturbative renormaliz- 
ability of the theory, and, in particular, confirms that the non-renormalizable perturbative 
counterterms play no special role in the asymptotic safety scenario. 

The rest of the paper is organized as follows. In Section [Til we review the counterterms 
arising from the perturbative quantization of general relativity coupled to a free scalar field, 
while in Section HTT1 we introduce the renormalization group methods employed. In Section HVl 
we revisit the Einstein-Hilbert truncation, for the pure gravity and matter-coupled cases, 
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and finally in Section |V] we present the results for our full fourth-order truncation. We 
conclude with a discussion of our results in Section [V]] All the details of the calculations are 
cointained in the three appendices: Appendix [A] contains the Hessians entering the FRGE 
for our truncation ansatz, Appendix Q3] presents the heat-kernel expansion for Lichnerowicz 
Laplacians, and finally Appendix O details the evaluation of the traces. 



II. PERTURB ATIVE NON-RENORMALIZ ABILITY AND COUNTERTERMS 



We start by reviewing the perturbative quantization of the Einstein-Hilbert action mini- 
mally coupled to a free scalar field. This provides the prototypical example of a gravitational 
theory which is perturbatively non-renormalizable at one-loop order as may be seen by 
computing its one-loop counterterms Ar dlv . In general, the one-loop effective action for a 
gauge theory takes the form 



r i-i°°P[$] = si$] + - STrln 



r2 ctot 

— (1) 



where $ A is the full set of fields (including auxiliary fields and ghosts), S ,tot [$] = S[<&\ + 
S* gf [$] + 5 ,gh [$] is the total action of the theory including the gauge-fixing and ghost terms 
S gi and S gh , and STr is a generalized functional trace carrying a minus sign for fermionic 
fields and a factor 2 for complex fields. Typically, this trace contains divergences which 
require regularization. 

Our starting point is the action 

%,</>] = J d A x^ K - 2 (_ J R + 2A) + ^^^0 , (2) 

supplemented by the gauge-fixing term 

S&i = ^J d ' x ^ r v F»F v , F, = D a h, a -ID^h, (3) 

and the corresponding ghost action. Here, k 2 = 167rG, G and A are the dimensionful 
Newton's and cosmological constant, respectively, denotes the Euclidean space-time 
metric, and is a real scalar field. The gauge-fixing is carried out via the background 
field method, splitting the metric and scalar fluctuations into a background part, g^, <f), 
and fluctuations around this background, /i^,/, according to g^ = g^ + h^ u and = 



4 



4> + f. Adapting the results obtained via the Schwinger-DeWitt technique, the one- 

loop divergences arising from ([2]) are readily found to be 1 

1 - [ d^xy/g \§R^ + iR 2 + Wo E + KW) 2 

' ' . (4) 

-K 2 {\R-2K)(d^<P) - f Ai? + 20A 2 , 

where e = (d-4) and £ = C^ p(T C^ ,7 -2i?^ J R^ + | J R 2 is the integrand of the Gauss-Bonnet 



Ar 



div 



(47r) 2 e 



term in four dimensions, with C^ up(T being the Weyl tensor. 

In order to get information on the renormalizability, the divergences 
considered on-shell. The equations of motion resulting from (j2J) are 



have to be 







R = 4A + \ k 2 (d^f 



R 



(5) 



Substituting these, eq. (TjJ can suggestively be written as 2 

1 



Ar 



div 



d x^fg 



213^ 203 . 463 ^ A 463 A „ 

E + R 2 RA + A 2 

360 80 20 10 



(6) 



As the R 2 and .E-terms are not of the form of the terms contained in the initial action, 
they cannot be absorbed by a renormalization of the coupling constants, indicating that 
the action (j2j) is indeed perturbatively non-renormalizable. The non-renormalizable on-shell 
counterterms are thus of fourth order in the gravitational sector and can be rewritten as 



Ar 



NR 



1 



(4vr) 2 e 



d^x^fg 



18 + 180 » vpa 



(7) 



There is a common prejudice (see, for example that these interactions have a devas- 
tating effect also on the possible non-perturbative renormalizability (asymptotic safety) of 
the theory Utilizing the new computational techniques developed in 28(, we will now show 
that this is not the case. 



1 There is a typo in the coefficient of the squared potential in [7[ , the correct formula is given in L ^ . 

2 Note that this expression agrees both with the one-loop counterterm found by 't Hooft and Veltman [l| 
for A = E = 0, and with the one of Christensen and Duff Sj once the contribution of the scalar field is 
subtracted. 
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III. THE FUNCTIONAL RENORMALIZATION GROUP EQUATION 



A powerful tool in the study of the renormalization properties of a theory is the Functional 



Renormalization Group Equation (FRGE) 31 1 

d t r fc [M] = |STr 



(8) 



where <3? denotes the physical fields and $ their background value. The FRGE describes the 
dependence of the effective average action $] on the coarse-graining (or renormalization 
group) scale k. Here, t = log(k/ko) and TZk{p 2 ) is a (matrix- valued) infrared cutoff which 
provides a fc-dependent mass-term for fluctuations with momenta p 2 < k 2 . Apart from the 
requirement that it interpolates monotonically between TZk{p 2 ) = as p 2 /k 2 — > oo and 
lZk(p 2 ) oc k 2 as p 2 jk 2 — *> 0, this cutoff can be arbitrarily chosen. For technical simplicity, 



37|, whose scalar part takes 



our subsequent analysis will be based on the optimized cutoff 
the form Rk(p 2 ) = ik 2 — p 2 )6(k 2 — p 2 ). 

The FRGE has two key features, owing mainly to the IR regulator structure. First, its 
solutions interpolate between the ordinary effective action V = and an initial action 

Ta at the UV cutoff scale, which in the limit A — > oo essentially reduces to the bare action 
(see |38] for more details). The effective average action is obtained by integrating out 
modes in the path integral from a UV cutoff scale A down to the scale k, as in a Wilsonian 
coarse graining procedure, with the modes below k being suppressed. Secondly, due to the 
derivative dtTZk{p 2 ) in the numerator, the contributions to the flow equation are localized 
on modes with momenta near k 2 , so that the trace remains finite and locally well-defined at 
all scales. In particular, while a theory might require a UV regulator at the level of its path 
integral, at the FRGE level this UV regularization is superfluous. 

The main shortcoming of the FRGE, however, is that it cannot be solved exactly. In order 
to extract physics from it, one therefore has to resort to approximations. One possibility 
is, of course, perturbation theory. In the one-loop approximation, where Tk under the STr 
is replaced by the fc-independent bare action, one then recovers upon integration the usual 



non-renormalizable logarithmic divergences 



27] 



Going beyond perturbation theory, a standard approximation scheme is the truncation of 
the RG flow, whereby the flow of the full theory is projected onto a subspace spanned by a 
finite number of interaction monomials. Making an ansatz r £.[<!>, = ^u^k) Oi[<&, $] with 



a finite subset of the interaction monomials 0% and substituting this ansatz into the FRGE, 
this technique allows one to extract the /3-functions for the dimensionful coupling constants 
Ui. When analyzing the properties of the RG flow, it is then most convenient to switch to the 
dimensionless coupling constants gi = k~ di Ui, with dj being the mass- dimension of Ui, which 
results in autonomous /3-functions d t gi = j3 gi (gi). Within Wilson's modern perspective on 
renormalization, the renormalizability of the theory is then determined by the fixed points 
(FP) {g*} of the /3-functions, {(3 gi (g*) = 0}. Around any such FP, the linearized RG flow 
d t gi = T$ij(gj — gj) is governed by the stability matrix 

By = <%/% } ■ (9) 

Defining the stability coefficients 9i as minus the eigenvalues of B, the relevant (irrelevant) 
directions are associated to the eigenvectors corresponding to stability coefficients with a 
positive (negative) real part. 

In general, it is useful to cast the effective average action into the form 

r fc [d>, $] = f fc [$] + f fe [$ -$,$] + S^[g,g] + S sh [g,g, ghosts] . (10) 

In this decomposition l\[$] depends on the physical fields only, and 5* gf and S sh denote 
the classical gauge-fixing and ghost-terms respectively encodes the deviations $ — $, 
thus vanishing for $ = $, and captures the quantum corrections to the gauge-fixing and 
ghost sector of the effective average action. For the remainder of this work, we will focus on 
truncations of the form 

f fe [$] = rr v M + r mattcr [^,0]. (ii) 

Here T^ rav is the gravitational part of the effective average action, for which we will specify 
two different truncations in Sec. [IV] and Sec. IVl while 

r matte %,0] = i J tfxjgg^d^d^ (12) 

is the (/c-independent) action for a minimally coupled free scalar field. Furthermore, we set 
Tfc [$ — $,$]= in the sequel. 3 For S fgf , we consider the following generalization of ([3]), which 



3 Of course, it would also be desireable to obtain a better understanding of the influence of — <£>, $] on 
the RG flow. In this context, the adaptation of the background independent version of Yk, discussed for 



Yang-Mills theories in 39(, could provide a valuable tool. 
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allows for a straightforward application to gravitational actions including higher-derivative 
terms 



cgf 

b 2 



l - J d'x^g F, Y^F U , F, = D v h^ - l -^D,h , Y*> = + pD 2 ]g^ . (13) 



The gauge-fixing (|2D is obtained as the limit p = 1, a = 1, (3 — 0. When analyzing the RG 
flows in sections [IV] and |V] however, it will be more convenient to set 

p = 0, a — > oo, (3 = 0, for the Einstein-Hilbert truncation, (14) 
p = 0, a = 0, (3 — > oo, for the higher- derivative truncation. (15) 

Here, the arrow indicates that the limit is to be taken under the trace of the flow equation. 

The key step in utilizing (jSj) for extracting /3-functions is the evaluation of the operator 
trace appearing on its r.h.s. For the ansatz (fTTj) . this STr decomposes into a trace over the 
gravitational and the matter sector, respectively. 



The contribution from the gravitational sector is obtained as follows 28] (for further 
details, see the Appendices). We first compute the second variation of T| rav and perform 
a transverse-traceless decomposition of the metric fluctuations and the ghost fields, dealing 
with the Jacobians arising from this decomposition by introducing suitable auxiliary fields 
as in the Faddeev-Popov trick. The subsequent computations can then be simplified by 
identifying g with a suitable class of background metrics g. This class must be general 
enough to distinguish the interaction monomials contained in r| rav and, at the same time, 
simple enough to ease the evaluation of the operator traces. For our present purposes, 
it suffices to consider the class of generic compact Einstein backgrounds (without Killing 
or conformal Killing vectors), not necessarily solving the equations of motion (J5]). Utilizing 
these backgrounds, all differential operators within our particular traces organize themselves 
into Lichnerowicz form 

A 2L </v = -D 2 V - 2R^<j> aP , A 1L M = [-D 2 - \R] cf>, , A OL ee -D 2 <f) , (16) 

i.e., minimal second order differential operators A s ^ = — D 2 + Q s , with spin-dependent 
matrix potentials Q s acting on transverse-traceless matrices (s = 2), transverse vectors 
(s = 1) and scalars (s = 0). This feature is crucial for the non-perturbative evaluation of 
the traces, as it makes them amenable to standard heat kernel techniques without having 
to resort to non-minimal (or fc-dependent) differential operators. The final steps in this 
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computation then follow the standard FRGE procedure (see e.g. [13|). First, the cutoff 
operators TZ s ,k are constructed in such a way that the modified propagators are obtained by 
replacing A sL — > -P s ,fc(A sZ/ ) = A sL + lZ St k(A sL ). Then, the traces in the flow equation are 
evaluated using the "early-time expansion" of the heat kernel adapted to the Lichnerowicz 
Laplacians. 

Written in terms of operator traces, the flow equation then takes the following generic 
form 

d t I\. [$,$] = S 2T + S hh + 5i T + S + n s 5 matter , (17) 

where n s gives the number of matter fields and the subscripts "2T", "IT" and "0" indicate 
traces taken on the space of symmetric transverse-traceless matrices, transverse vectors and 
scalars, respectively. Applying the background-field method (setting = 0, for convenience), 
it is straightforward to find the matter contribution to the flow equation 32J , 



'matter 



:Tr 



d t R 



P, 



(18) 



Furthermore, owed to the special gauge choices ffl"3l) with ffT4|) and ffT5l) . Sn and So are 
universal, in the sense that they are independent of the particular rf rav [g] adopted here, 



S 



IT 



— Trix 
2 1T 



Pi,k 



So = -2 Tr ° 



d t R 



3-Po,fc — R 



(19) 



Following the computations outlined in Appendix [Cj the evaluation of the traces can 
be carried out using the early-time heat-kernel expansion for Lichnerowicz-operators. The 
resulting expressions are given in eqs. (1C9j) and (1C10I) . respectively. We should also note 
that it is straightforward to consider the effect of a number n s of scalar fields by adding 
further copies of (1C9I) to the RG equations. In what follows, however, we will mostly focus 
on the cases n s = 0, 1, keeping the label n s only to highlight the matter contribution. 



IV. FIXED POINTS OF THE EINSTEIN-HILBERT TRUNCATION 



In this section, we approximate r| rav by the Einstein-Hilbert action with scale dependent 
coupling constants, 

= TTrV / d*Xy/g{-R + 2A k ) . (20) 
lo7rGfc J 

This truncation has already been considered in the context of pure gravity in 17|, [l8|, [l9|, |20(] 



pgrav p 
1 k 



and in the matter-coupled case in 



27|, 



32|, 



351 ] . Here, we complement this analysis 



by implementing the gauge-fixing (I13H14I) and organizing the operators inside the trace in 
terms of Lichnerowicz Laplacians on a generic Einstein space, lending further evidence to 
the robustness of these earlier works. 

Using the results of Appendix |A], the non-universal traces resulting from f[2"Uj) are 



<S 2 T — -Tr 2 T 



d t (uiR k 



uiP 2 ,k + uxR/2 + m 



Shh — oTro 



d t (uiR k 



2 [3uiPo !k + 2u Q _ 



(21) 



where the coupling constants Ui are defined in ( l27|) . Adapting the computations outlined in 
0, 122I to the steps described in Section lITlj the /^-functions of the dimensionless Newton's 
constant g k = G k k 2 and cosmological constant X k = A k k~ 2 then read 

(3 g =(2 + , (3 X = (qjf - 2)A + |-S 3 (A) - J^ H (s^(-2A) + S*(-f A)) , (22) 

where, 

eh _ 9Bi{\) , . 

Vn ~ Gtt - gB 2 (X) (23) 

denotes the anomalous dimension of Newton's constant and where we have defined 
B 1 (X) = K - 7)$1(0) - 2^(0) + - 10$1(-2A) - 15$|(-2A) , 

B 2 (X) = 5*J(-2A) + f m-2X) - |^(-|A) , (24) 
B 3 (X) = K - 4)^(0) + ^(-|A) + 5^(-2A) . 

We can see that the inclusion of matter fields simply results in a shift by a constant in 
the equations above. In a sense, "small" values of n s could therefore be interpreted as a 
"perturbation" of the /3-functions for pure gravity, which are recovered in the limit n s = 0. 
Note that these /3-functions are non-perturbative, in the sense that the anomalous dimension 
7/™ contains infinitely many powers of g. 

Analyzing the fixed point structure of the /3-functions ( T22|) . we first note that both the 
one-loop and the non-perturbative /3-functions possess a GFP at A* = 0, g* = for all values 
of n s , corresponding to the free theory. Its stability coefficients are given by the canonical 
mass-dimensions of A and G, that is, Q\ — 2 and 9 2 = —2. Owing to the negative mass 
dimension of Newton's constant, the only relevant direction is at G = 0, amounting to a 
trivial theory with just the cosmological term. As soon as we turn on Newton's constant, 
the flow is carried away from the GFP in the UV. Thus, our gravity-matter theory is not 
inside the UV critical surface of the GFP, verifying its perturbative non-renormalizability 
from the Wilsonian viewpoint. Note that this behavior is independent of n s . 
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Ref. 


9* 


A* 


g*X* 


& ± iO" 


a 


p 


cutoff 


Here 


0.902 


0.109 


0.099 


2.52±1.78i 


CO 





II, opt 


RS 


0.403 


0.330 


0.133 


1.94 ± 3.15i 


1 


1 


I, sharp 


LR 


0.272 


0.348 


0.095 


1.55 ± 3.84i 


1 


1 


I, exp 




0.344 


0.339 


0.117 


1.86 ± 4.08i 


CO 


1 


I, exp 


L 


3vr/8 


0.25 


3vr/32 


1.67 ± 4.31i 


CO 


1 


I, opt 


CPR 


0.707 


0.193 


0.137 


1.48 ± 3.04i 


1 


1 


I, opt 




0.556 


0.092 


0.051 


2.43 ± 1.27i 


1 


1 


II, opt 




0.332 


0.274 


0.091 


1.75 ± 2.07i 


1 


1 


III, opt 



TABLE I: Einstein-Hilbert truncation: comparison of the characteristics of the pure gravity NGFP 
obtained here employing the "universal gauge-fixing" , and results reported earlier in the literature 

H 



IJWLR), 



|2fl|| (L), and [27j (CPR), respectively. The cutoff classification as Type I, II 
or III follows [27], while with "opt", "sharp" and "exp" we refer to the shape function being of the 
optimized, sharp or exponential type respectively. The fixed point is robust under variation of the 
gauge-fixing parameter a, the shape function used in the IR regulator, and the implementation of 
the regularization. 



Remarkably, the /3-functions fl22|) also give rise to a NGFP at positive values A* > 0, g* > 
0. Its position and stability coefficients for n s = 0,1 are given in Tables HI and ILT1 (together 
with a comparison to earlier works). Note that this NGFP is UV-attractive for both the 
dimensionless Newton's constant and the cosmological constant. The gravity- matter theory 
considered here is within the UV critical surface of the NGFP. In other words, matter- 
coupled gravity in the Einstein-Hilbert truncation is asymptotically safe. Note also that 
making the transition from pure gravity to gravity coupled to a free scalar field has a rather 
small effect on the numerical values obtained for the NGFP. 



V. FIXED POINTS OF THE HIGHER-DERIVATIVE-MATTER TRUNCATION 

The key question raised by the results obtained within the Einstein-Hilbert truncation is 
whether the resulting fixed points survive in the full theory and, in particular, whether the 
NGFP will persist, with similar characteristics, once the perturbative counterterms ([7]) are 
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Ref. 


9* 


A* 


9* A* 


9' ± W" 


a 


p 


cutoff 


Here 


0.860 


0.131 


0.112 


2.58 ± 1.95i 


oo 





II, opt 


PP 


0.254 


0.366 


0.093 


1.71 ±4.16i 


1 


1 


I, exp 




0.320 


0.359 


0.115 


2.08 ± 4.38i 


oo 


1 


I, exp 



TABLE II: Einstein-Hilbert truncation: comparison of the characteristics of the NGFP obtained 
from gravity coupled minimally to a real scalar field. The first line is obtained from the "universal 
gaugefixing" , while the data of the second and third line has been obtained in 



34, 



35|] (PP), and 



is given for further comparative purposes. Again, the fixed point is robust under variation of the 
gauge-fixing parameter a, the shape function used in the IR regulator, and the implementation of 
the regularization. 



included in the truncation subspace. While the former is a million- dollar question (recession 
notwithstanding), the latter can be answered positively. 

To this effect, we enhance the truncation subspace (J2"0j) and consider the ansatz 



-pgrav 
1 k 



[g] 



-R + 2A A 



1 



0i 



R z + —C^ pa C^ + — E 



(25) 



167rG fc v " 7 3a k ' 2a k ■ ■ a k 

which is precisely of the form Einstein-Hilbert-action plus perturbative counterterms. In the 
spirit of the RG, the numerical coefficients in the latter have been replaced by the canonical 
(scale-dependent) coupling constants. 

Following the derivation given in Appendix |2J the gravitational contribution to ifPTj) 
enters into S 2 t and Shh only and is given by 



S 



2T 



~>hh 



1 

— Tr9T 

2 



2 



d t {2u 3 (P, 



2,k 



(Ml 



+ u\,R)R 2 ,k} 



2u 3 Pl k 



ui + u\>R)P 2j k - \u\R - uq 
8 t {6« 2 (P % - A 2 0L ) + ( Ul - 2u 2 R)R 0yk } 



(26) 



Qu 2 p o,k + Oi - 2u 2 R)P ,k + f MO 
The coupling constants appearing in these expressions are related to (|25l) via 
A k 1 



Uq 



&nG k 1 \QnG k z 3a k Qa k 7 

and u\) = 2u 2 — Note that, because of the Einstein-space choice, we can distinguish 
only two of the three higher- derivative couplings. Lastly, note also that including 5 matt er 



u 2 



<~^k 6k 



U 3 



1 

2a k o k 



(27) 



12 



has a similar effect as in the Einstein-Hilbert case, leading to shifts in certain coefficients 
appearing in the /3-functions. 

The projection of the traces onto the truncation subspace spanned by can again be 
carried out utilizing the early time expansion of the heat kernel adapted to the Lichnerowicz 
operators on a general Einstein background, as detailed in Appendix [B] Following the 
computation outlined in Appendix O, and introducing the dimensionless couplings 

g Q = u k~ 4 , gi = u x k' 2 , g 2 = u 2 , # 3 = «3 , (28) 

the /^-functions following from our truncation are given by 

d t go = - 4£o + {Ci + C, + (2n s - 8)<^} , 

1 s~ „■„ (29) 

dm = (4^)2 { sio^ 3 + i^ 3 + 12 ^'^} ' 

where the threshold functions ip are respectively defined in (1C3P and (107j) . with the for- 
mer evaluated at zero argument. The expansion coefficients Cj and Cj arise from evaluating 
the traces S 2 t and Shh, respectively, and are defined in eqs. ( 1C14I) and ( 1C16I) . For notational 
reasons, the /3-functions (1291) are given implicit form. In particular, we stress that both the 
left and right-and-side contains derivatives d t gi- The "standard" /3-functions d t gt = Pi(gi) 
can then be obtained by solving these equations for dtgi, which can be straightforwardly 
done using algebraic manipulation software. 

The resulting expressions may again be expanded for small g and, here, a. In this respect, 
we first note that ( [291 contains contributions from arbitrary powers in g, a, and hence that 
the /^-functions capture some truly non-perturbative information. Secondly, we have verified 
that the leading contributions in this expansion reproduce the known universal parts of the 
one-loop /3-functions in higher-derivative gravity, providing an important confirmation of 
the correctness of our derivation. 

Remarkably, the fixed point structure originating from these higher-derivative /3-functions 
is very similar to the Einstein-Hilbert case. First, we recover the two generalizations of the 
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GFP, familiar from perturbation theory 4 , 

g* = 0, A* = 0, ct* = 0, u>1 )2 = -± (90 ± V15708 - 101n s ) , (30) 
existing for < n s < 155 and with stability properties given by the following eigensystem 

01 = 2, K = {1,0,0,0} T , 6 2 = -2, K 2 = {^,1,0,0} T , 
03 = 0, U 3 = {0,0,1,0} T , 04 = 0, V4 = {0,0,0,0} T . 

These GFPs correspond to the free theory, and their stability coefficients are given by the 
canonical mass dimension of the corresponding (dimensionful) couplings. In particular, the 
eigendirection associated with Newton's constant is still UV repulsive, while the directions 
associated with the new couplings <Jk,^k are marginal. Going beyond the linear approxima- 
tion, the marginal directions are found to be UV-attractive, in accordance with the one-loop 
calculations [401 ] . 

Most importantly, the matter-coupled higher-derivative truncation also gives rise to the 
generalization of the NGFP. Its corresponding position and stability coefficients are given 
in Tables HTT1 and UVl (under the entries "i? 2 + C 2 +scalar"). For completeness, these tables 
also include the data on the NGFP for the pure gravity case ("i? 2 + C 2 "), first reported in 



281 ] . and we note that its properties are again very similar to those in the gravity- matter 



case, thus giving rise to essentially the same picture. 



stability co- 



40j, but it is 



One salient difference with the Einstein-Hilbert case is the fact that al 
efficients are now real. This is in agreement with the one-loop results of 
surprising that, unlike in the Einstein-Hilbert case, the transition from the one- loop to the 
non-perturbative treatment does not give rise to complex eigenvalues. We can trace this 
result to the contribution of the C 2 terms coming out of the traces: indeed, restricting our 
computation to a spherically symmetric space we again find complex eigenvalues. 

Crucially, increasing the dimension of the truncation subspace, with respect to the 
Einstein-Hilbert case, adds one UV-attractive and one UV-repulsive eigendirection to the 
stability matrix, so that the UV critical hypersurface in the truncation subspace is now 
three-dimensional. We then have a three-dimensional subspace of RG trajectories which 



4 In a slight abuse of notation, we neglect the topological term here,setting 6 = 0, dtO — 0. Also note that 
the existence of this GFP is actually compatible with the analysis [24j , which did not consider the case of 
the inverse coupling 1/a — > oo. 
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Truncation 


g* 


A* 


U2 


«3 


U4 


X*g* 


cutoff 


R 2 + C 2 


1.960 


0.218 


0.008 


-0.0050 




0.427 


II, opt 


LR II 


0.292 


0.330 


0.005 






0.096 


I, exp 


CP 


1.389 


0.221 


* 


* 


* 


0.307 


perturbative one-loop 


R 2 + C 2 + scalar 


2.279 


0.251 


0.010 


-0.0043 




0.571 


II, opt 



TABLE III: Position of the NGFP optained from the non-perturbative /3-functions of the R 2 + C 2 
truncation, eq. f|29|) . For comparison, we also give the data of the i? 2 -truncation 
the perturbative one-loop result 



2J (LR II), and 
= -0.0228, 9* = 



40J (CP). In the latter line, the * indicates that lo* - 
0.327 approach finite values in the UV, while a runs logarithmically to zero, realizing the asymptotic 
freedom of the one-loop result. The last line gives the position of the NGFP upon including a 
minimally coupled scalar field. 



Truncation 


01 


02 


03 


0A 


R 2 + C 2 


2.51 


1.69 


8.40 


-2.11 


LR II 


2.15 + 3.79i 


2.15 - 3.79i 


28.8 




CP 


4 


2 


* 


* 


R 2 + C 2 + scalar 


2.67 


1.39 


7.86 


-1.50 



TABLE IV: Stability coefficients of the NGFP optained from the non-perturbative /3-functions of 
the R 2 + C 2 -truncation, eq. (|29p . For comparison, we also give the data of the i? 2 -truncation 



40] (CP). In the latter line, the * indicates the 



24] (LR II), and the perturbative one- loop result 
logarithmic running of the marginal coupling constants towards asymptotic freedom. The last line 
gives the stability coefficients of the NGFP upon including a minimally coupled scalar field. 



are attracted to the NGFP in the UV and are therefore "asymptotically safe". Thus, non- 
perturbative renormalizability persists also in the presence of the one-loop perturbative 
counterterms in the truncation ansatz. 



15 



VI. DISCUSSION AND CONCLUSION 



In this paper, we have analyzed the fixed point structure underlying the renormalization 
group (RG) flow of gravity minimally coupled to a free scalar field, within a truncation 
approximation. From the viewpoint of perturbative quantization, this setup provides a pro- 
totypical example of a quantum theory of gravity which is perturbatively non-renormalizable 
at the one-loop level [lfl. Here, higher-derivative interactions arise as perturbative counter- 
terms, signaling the presence of divergences which cannot be absorbed by the renormalization 
of the coupling constants. However, despite the breakdown of the perturbative quantization 
scheme, there is the possibility that this gravity-scalar theory constitutes a well-defined and 
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35|. With this in 



predictive quantum theory within the realm of asymptotic safety 
mind, we first considered the case of the Einstein-Hilbert truncation, before extending it 
to a higher- derivative truncation by including the interactions of the form of the one-loop 
counterterms. 

As our main result, we show that all these truncations give rise to a non-Gaussian fixed 
point, which underlies the conjectured asymptotic safety of the theory, in addition to a 
Gaussian fixed point linked to the perturbative quantization. Both fixed points are robust 
under the extension of the truncation subspace by higher derivative terms. This result ex- 
plicitly shows that, contrary to a common worry, the inclusion of perturbative counterterms 
in the truncation subspace of a gravity-matter theory has no qualitative effect on its fixed 
point structure. In particular, we find no indication that these interactions are fatal to 
non-perturbative renormalizability of the theory. 

A second remarkable property surfaces when comparing the fixed point structure obtained 
for pure gravity (n s = 0) and gravity coupled to one free scalar (n s = 1) given in the top 
and bottom lines of Tables HTT1 and HVl respectively. Including the scalar field shifts the fixed 
point values obtained for pure gravity only very mildly, so that the resulting fixed point 
patterns are very similar. In a sense, this indicates that (at least for the present truncations) 



the UV behavior of the gravity-matter theory is still dominated by its gravitational 



sector, 



so that it still behaves "essentially gravitational" at high energies. Following [34], |4l|, it 
would be very interesting to determine which matter sectors lead to asymptotically safety 
gravity- matter theories (which we might dub the "asymptotic safety territories"), taking the 
higher- derivative terms (J25l) into account. 



16 



While our results on the interplay between the perturbative counterterms and asymp- 
totic safety in the gravity-matter case are already trend-setting, it would nevertheless be 
desirable to carry out an analogous computation for pure gravity, where non-renormalizable 
divergences set in at two- loop level 3|, |4|. This is, however, still beyond the current 
technical scope of the functional renormalization group techniques employed in this paper. 
Nevertheless, various arguments have been put forward 
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42] that the situation there will 



be similar to the one encountered here: perturbative counterterms are likely to have no 
special effect on the asymptotic safety of the theory. 
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APPENDIX A: THE HESSIAN FOR HIGHER-DERIVATIVE GRAVITY 

In four dimensions, the derivative expansion of r| rav [(7] up to fourth order can be organized 
into the five interaction monomials, 

I = J d 4 x^/g , h = J d A x^gR , 

(Al) 

h = J d 4 x^/gR 2 , I 3 = / d*x^gR pv R^ , h = J d A x^/gE 

where E = R 2 — AR llv R tJ,u + R^ upa R^ vpa is the integrand of the Euler topological invariant, 
J d A x^fg~E = 32tt 2 x- Constructing the argument of the traces entering into the FRGE 
requires the second variation of these invariants. In this context, we first note that I4 
is a topological quantity, so that its variation with respect to the metric vanishes. To 
obtain the hessians of the other invariants, we split g^ v = g^ v + h^ u , where g pu denotes 
a fixed background metric and is an arbitrary fluctuation. The genera 
'or these variations, valid for an arbitrary background g^ V) can be found in 



expressions 



46, 



47J, 



48|, 



451 ] (see also 



49]). For our purposes, however, it suffices to consider these variations on 
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backgrounds = g? v , where the index £ indicates that the background metric is a generic 
Einstein metric. These are metrics satisfying R^ u = 4<? Mi , (but not necessarily C^ vp(T = 0) 
and, using the contracted Bianchi identity, this condition also implies that D x R\ apu = 0. 
For these spaces, the Hessians of I n then simplify considerably. At the two-derivative level, 
we obtain 



5 2 I = j^x^[\h 2 -\h^] , 

5 2 h = j d A x^g~[\h a P [D 2 - \R] h af3 + R^h^W 

- \hD 2 h + h(D a D^h af3 ) + (D^ a )(D„h" 
while the variations of the four-derivative terms yield 



(A2) 



5 2 h 



j rf 4 xv^| 2h [D 4 - ±R 2 ] h + Rh aP D 2 h aP + 2RR a ^ v h aP W lv + 2{D a D j3 h 
+ (D a D p h a P) [-AD 2 + R]h + 2R(D a h^)(D"h^)} , 



a/3\2 



(A3) 



and 



5 2 / 3 = J d A x Vg~{ \h a ? [D A + \RD 2 ] h a p + \h [D A - \RD 2 - \R 2 ] h 



D 2 -% 



h + (D a h 



D 2 + 3_R 



vP, 



(A4) 



+ {D a Dph af) ) 2 - 2h aP R a ^ [D 2 + \R] h, u + 2h a pR aX ^R^ au h^ 



respectively. Here, the bar denotes that the corresponding quantity is constructed from the 
background metric and h = g^h^. 

A remarkable feature of these variations is that they can naturally be written in terms of 
second order minimal operators of Lichnerowicz form (I16p . In particular, the four- derivative 
operators appearing in ( 1A4I) and ( 1A3I) factorize into squares of these (modified) Laplacians. 
Performing the TT-decomposition fIBlOl) for the metric, a brief computation establishes 



S 2 I 
S 2 h 



J d A x^{lh 2 - \h^hl u - £"A ir £„ - \a [3A 0L - R] A 0L a} 



d 4 xV§\MA 0L h-lh T ^ 



A 2L + \R 



(A5) 



+ \h [3A 0L - R] A 0L a + ±a [A 0L ~ R] [3A ol - R] A 0L a} 
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For the four-derivative terms, an analogous computation shows 

5 2 I 2 = J | h [3A 0L -R]A 0L h-R h a ^ A 2L h r ap 

+ fa [3A oi - R] Aq L a + lh [3A 0L - R] A 2 0L a} 



(A6) 



and 



5 2 h = J £ d 4 x^{ \ h Ta ? [A 2L - \R] A 2L hl p + \h [3A 0L - R] A 0L h 
+ ±<r [3A 0L - R] A 3 0L a + \h [3A 0L - R] A 2 QL a} . 



(A7) 



As a welcome side-effect, we also observe that the introduction of the Lichnerowicz- 
Laplacians diagonalizes the transverse-traceless sector of the fluctuations. With respect 
to "off-diagonal" terms in the metric sector of r? [y], the R 2 + C 2 -truncation has thus 
the same level of complexity as previous computations which included (polynomials of) the 
Ricci-scalar only and referred to a maximally symmetric background. 

In order to complete the construction of the operator traces, we now turn to the gauge- 
fixing and ghost terms originating from ( JT3l) . For the higher-derivative action of Section IVl 
we thereby work with a = 0, p = 0. In this case, the TT-decomposition of S gf yields 

S gf = -f / ^v^{e4( A ^ + f) A ye + ^[(!AoL-f) 2 (AoL-f)A OL ].7} (A8) 

The ghost sector now contains, in addition to the usual (complex) C, C-ghost fields, a third 
ghost 45| due to the two-derivative contribution (det (3D 2 ) 1 / 2 . Introducing the complex- 
valued Grassmann fields B^, B^ and the real field for the latter term, and TT-decomposing 
the ghost sector of the resulting action then leads to 

Sc-thost = - J d 4 xVs{c^A 1L C^ + \n[^oL - R] Aolv} , 

S£SL = - / [Ail + f ] B T " + B [A 0L - f ] A 0L B (A9) 

+ ¥l[^L + W» + \b [A 0L -|]A 0L 6}. 

Note that, in the literature on higher-derivative gravity, the contribution of the B, 5-ghost 
field is usually absorbed into the usual C, C-ghost, hence the need of only a third (real) 
ghost. We prefer here to introduce a fourth ghost to clearly separate the higher-derivative 
contribution from the usual second order term. The two choices are of course equivalent. 
In the following, we impose a "mode by mode" cancellation between the gauge-degrees of 
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freedom in the metric and the ghost sector [25( , which results in a precise cancellation of all 
the "unphysical mode contributions" to (fl9j) . 

Finally, there are additional contributions to the flow equation arising from the Jacobi- 
determinants introduced via the TT-decomposition, 

1 /2 

Jgrav = (det' (1Ti0) [M^] ) , Jc- g h = J B - gh = (det'^])- 1 , J b _ gh = (det'[A 0L ]) 1/2 . 

(A10) 

Here, the primes indicate that the unphysical modes are left out from the determinants. 
Furthermore, M^'^ is a (d + 1) x (d + l)-matrix differential operator whose first d columns 
act on the transverse spin one fields £^ and whose last column acts on the spin zero fields a 
and which reads 



R r>v 3 A 2 R A 
2 U 4 A 0L ~ J A 0L 



(All) 



In order to account for these contributions, we follow earlier works [25|, |26|, |27| and introduce 
appropriate auxiliary fields so as to exponentiate these determinants via the Faddeev- Popov 
trick. The resulting "auxiliary action" then becomes 

ST*= [ ^^{[Cj^][M^]'[Cj,a;] T + [ C -J,c][M^)]lc^c] T 

J 1 (A12) 

+ s A' 0L s + tA' 0L t + x Kl X + \<t> Aol 0} • 

Here the gravitational sector contains the transverse ghost cE, c Tm , a "longitudinal" Grass- 
mann scalar c, c, a transverse vector £j and a real scalar 00, while the ghost determinants are 
captured by the contribution of the complex scalar fields s, s, t, i, the complex Grassmann 
fields Xi X) an d the real scalar field (p. 

(2) 

We now have all the ingredients for constructing all the Hessians T k required in the 
r.h.s. of the flow equation. These are collected in Table |V] The gravitational sector then 
arises from combining the contributions from the variations (1A5jl . (1A6I) . and ( 1A7I) with the 
gauge-fixing action ( 1A8j) . reinstalling the corresponding coupling constants. These kernels 
contain all the information required for constructing the r.h.s. of the flow equation (fT7|) . To 
cast the result into the form (|26|) and (fT9|) it is thereby useful to note that the <r-/i-crossterm 
vanishes in the limit (3 — > 00. Thus, the combined contribution from a and h splits into the 
sum of the /i/i-trace (126]) and the contribution of the acr-part. The latter can be combined 
with the contribution of all the other scalar fields to give rise to the universal scalar trace 
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Index 


(2) 

Hessian T y k ' 




2u 3 A| L - [2u 2 R - lu 3 R + ui]A 2 l - \u\R - u 




-/3[A 1L + f ] A{ L - [ Ul R + 2u ]A 1L 


7 7 

hh 


1 l~i o A 2 i o / c\ I~> \ A i o i 

^ [18u 2 A5 L + 3(tii - 2u 2 R)A 0L + 2uq\ 


era 


-fg [3A l - R] [A 0L - f ] Aol + | [6u 2 A§ L + ^Aol - - 2u ] [3A 0L - A l 


ha 


| [ Ul + 6 u 2 A 0L ] [3A 0L - -R] A 0L 




/3Ai L 


fjrj 


f [3 A 0L Aol 


Bi B 1 » 


-/3[Ail + |] 


BB 


-/3[A l - f ]A l 


iTiTii 

bf. b 1 v 


-/3[Ail + f ] 


bb 


-/3[A 0L - f ]A 0L 




4A 1L 


UJU) 


1 rt-> A 7~»1 A 

i [3 Aol - R] A l 




2A 1L 




IfSAnr — RlAnr 


00 


Aol 
Aql 



TABLE V: Matrix entries of the operator 1^ in the gravitational, ghost and auxiliary sector 
(separated by the horizontal lines), respectively. The elements are symmetric under the change of 
bosonic indices, while they acquire a minus sign when Grassmann-valued indices are swapped. 

Sq. Similarly, combining all contributions from transverse vectors leads to iSit, which is also 
independent of the details of the gravitational action. Finally, the /i T /i T -sector produces the 
52T-trace. 

Lastly, we note that the derivation of the flow equation for the Einstein-Hilbert truncation 
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proceeds in an entirely analogous manner. In this case, the gravitational sector arises from 
the contributions of flA5j) with the TT-decomposed gauge-fixing term f|T3|) in the limit a — ► 
oo, (3 = and p = 0, with a similar vanishing of the <7-/i-crossterm and decoupling of the 
h and a traces. The ghost sector now contains only the C-ghosts, and the auxiliary sector 
consequently does not contain the 0, x an d t fields. Combining the aa with all the other 
scalar field traces and the with all the other transverse vector traces then results in the 
universal traces Sq and <Sit, respectively, which are again given by ffT9l) . The remaining 
S 2 t and Shh traces can then be straightforwardly constructed from Table [V] by setting the 
higher-derivative couplings to zero. This concludes our derivation of the gravitational sector 
of the flow equation (JTTj). 



APPENDIX B: HEAT-KERNEL COEFFICIENTS FOR LICHNEROWICZ 
LAPLACIANS 



For evaluating the operator traces appearing in Section IV] we require the heat-kernel 
expansion for the Lichnerowicz operators (TTB1) . evaluated at a generic Einstein manifold, up 
to fourth order in the derivative expansion. In this appendix we derive the corresponding 
coefficients starting from the early time heat-kernel expansion of a generic two-derivative 
differential operator 43|, |44| (see also 27j for a nice exposition in the context of the FRG). 



1. Heat-kernel coefficients for unconstrained fields 



In general, the early time heat-kernel expansion of a generic second order differential 
operator A = —D 2 + Q takes the form 



Tr [e itA ] = (r^J J d 4 x^g {tr a - it tr a 2 - t 2 tr a 4 + • • •} 



(Bl) 



with the heat-kernel coefficients a^k given by 43 1 

a = 1 , a 2 = P , 

1 111 ( B2 ) 

a 4 = — (R P vapR^ ~ R»uR» v + D 2 R) 1 + -P 2 + —Kyv-RT + -D 2 P . 

loll Z LZ 

Here, D 2 is the covariant Laplacian with respect to the (background) metric, Q is a matrix- 
valued potential, P = ^Rl + Q, IZ^ = 2D^ t D u ] is the commutator of the covariant deriva- 
tives, and tr denotes a trace with respect to the spin-indices of the fields on which A acts. 
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For the purpose of this paper, we have to evaluate tr s a%k for scalars (s = 0), vectors (s = 1), 
and symmetric tensors (s = 2). In the latter two cases, the tr s are defined as 

respectively. The matrices r R, lLV W ,v are trivial for the scalar case, whereas for vectors and 
tensors they respectively read 

[n a/3 n af3 ] = - R a p lfl R a ^ u , 

(B4) 

-i-> fp a ^~\ — p pa/37 n p p a /37 „ 1 op P? a <^ 

The differential operators appearing in the traces (J26l) and (fl9l) are the Lichnerowicz 
operators f[T6"j) . i.e., second order differential operators with matrix-potentials 

Qo = , [Ql]/iiv = -^Q^uR-, [Q,2\nvaP — ZRfiavP • (B5) 

Their heat-kernel coefficients on a generic four-dimensional Einstein manifold without 
boundary can be obtained by substituting these potentials into the expressions for the generic 
heat-kernel expansion. Evaluating the spin-traces, we obtain 

troflo = 1 j troa2 = qR , tr ci4 = ■^QRp Ua f3R fJ ' ual3 + ^-^ 2 > 

triao = 4, tria2 = |-R, tr^ = — -^RpvapR^^ + j^R 2 > (B6) 

ti2 a o = 10 , tr2a2 = \R 1 tr2a 4 = j^Rpv a pR^ ua ^ ~ ~^iR 2 ■ 
This result completes the heat-kernel expansion for unconstrained fields. 



2. Heat-kernel coefficients for fields with differential constraints 

In order to apply the early-time heat-kernel expansion to the operator traces (126]) . and 
( fj~9l the heat-kernel coefficients for the unconstrained fields given in the last subsection must 
be converted into the expansion coefficients for the transverse vectors (IT) and transverse- 
traceless symmetric matrices (2T) entering into the TT-decomposition. 

In the decomposition of a vector field into its transverse and longitudinal parts, 

A, = Al + Dp$, D"AJ = 0, (B7) 
the spectra of -D M $ and $ are related by 

A 1L D^ = Dp(A 0L - , (B8) 
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and the constant mode in $ does not contribute to A^. Thus, the decomposition of the 
s = 1-trace takes the form 

Tr 1T [e itA ^] = Tn [e itAli ] - Tr [S A ° L ~ R ^] + e~ itR / 2 . (B9) 

where the last term removes the constant $-mode from the s = 0-trace. A similar argument 
applies to the TT-decomposition of the symmetric tensor, 

1 1 

V = V + 2D (^) + D ^ D v a + ^Qfw a olct + -g^h , (BIO) 

where the components appearing on the RHS of this decomposition are subject to the 
constraints 

<T/£ = 0, D»h^ = 0, D^ = 0, h = g, u h^. (Bll) 
In this case, one can use 

A 2L [D„D V + \g^A QL ] a = [£>„!>„ + J^Aql] [A ol - f ] a , (B12) 
A2L9itvh = Qfiu [Aql - f ] /i , 
to relate the spectrum of A 2 ^ to the ones of the vector and scalar fields. Furthermore, (IBlOj) 



indicates that the constant mode in a, scalars subject to [D fl D u + j^Aoz,] a = 0, and 
transverse vectors satisfying = do not contribute to h^ u , so that the corresponding 

modes have to be removed from the decomposed spectrum. By contracting the last two 
equations with D", one can show that these are eigenmodes of A l and Ail with eigenvalues 
Aol = 0, A l = -|, and A 1L = 0, respectively. 5 The multiplicity of the latter two is given by 
the number of Killing vectors %v and conformal Killing vectors tickv of the background. 
Taking into account ( 1B9|) . the operator trace for transverse-traceless tensors field can then 
be expressed in terms of traces over unconstrained fields 

Tr 2T [e" A -] = Tr 2 [e itA -] - Tr x [e itA -] - Tr [S A ^ R ^] + n KV + ncxve"^ 6 . (B13) 

In the following, we will assume that our background is generic, in the sense that its metric 
does not admit Killing or conformal Killing vectors. 



For a spherical background, these coincide with the two lowest eigenmodes of —D 2 acting on scalars and 



the lowest eigenmode of — D 2 acting on vector fields [24 \27\ 
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From eqs. (1B9j) and (1B13I) it is then straightforward to compute the heat-kernel coefficients 
for Lichnerowicz Laplacians acting on transverse vectors and transverse traceless symmetric 
matrices. For a generic Einstein background, these read 



These coefficients are the key ingredient for evaluating the operator traces (1261) and (fl9|) and 
constitute the main result of this appendix. 

APPENDIX C: OPERATOR TRACES AND /^-FUNCTIONS 

In this appendix, we evaluate the operator traces appearing on the r.h.s. of eq. (fTTI) . We 
start with reviewing some general properties and definitions before computing the traces 
entering into our truncations explicitly. 

1. General trace technology 

The key observation for evaluating the traces entering f[T7|) is that they contain only 
minimal second order differential operators, which commute with all other elements (like 
the curvature scalars) inside the trace. Their projection onto the truncation subspace can 
then be found using the heat-kernel coefficients for constrained fields given in (IB14I) . Here, 
the key formula is 



Tr[W(A lL )] = — — / d 4 x^{Q 2 [W]tx i a + Q 1 [W}tr i a 2 + Q [W}tr i a± + ...} , (CI) 



where W(z) is a smooth function whose argument has been replaced by the Lichnerowicz 
operators and where the dots indicate higher-derivative terms at order six and higher, which 
are outside our truncation subspace. The functionals Q n [PF],ra > are defined as 



tro^o — 1 , tr a2 — \R 



, tlQ(l4 — R^ ua pR^ @ + gQ R , 

, tr^G^ = — Rp Va pR^ va ^ + <jfj) R 2 j 

5 p „ 10 p TDauaS 29 p2 



tr^'p'^o — 3 , tr^ r pd2 — R , 



(B14) 



tf2T a — 5 , tr2T^2 — ~o 






(C2) 
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In order to construct the /3-functions for the dimensionless couplings, it is useful to convert 
the QnfPF] into standardized dimensionless threshold functions, 

1 



r(n) Jo 
1 



n{ ] r n) J 



dz z 



(z + RW(z)+w)p 
RW(z) 



(C3) 



(z + R(°)(z) +w)p ' 
and their generalizations for higher-derivative theories, 

i r°. n _! (z+ J R(°)(^)) m (i?(°)(z)-^ o )'(z)) 



T'Ju, v, w) := 



dz z 



1 



dz z 



u(z + R(V(z)y + v(z + R(°)(z)) + w 
n _ x (z + R^(z)) m (2z + R®(z)) 1 RW(z) 



p ■ 



p ■ 



(C4) 



r H ( z + R(o)( z )) z + v{z + R(P)( Z )) + w'^ 

defined for n > 0. For the particular functions W occurring in this paper, these relations 



are 



Q, 



9t(gkRk) 
{2g k ){P k +c k )P 



k 2(n- P+ i) fofa/tf) + l dt l n (^) <S> p (c k /k 2 )) , n > 0, 



(C5) 



in the case of the Einstein-Hilbert truncation, while for the R 2 + C 2 -truncation we addition- 
ally have 



Q, 



d t (g k (P 2 -A 2 ) + ~g k R k 
(u k P 2 + v k P k + W k )P 



= n ~ 2p+2 ^ {d t g k f p nfiA (u k ,v k /k 2 ,w k /k 4 ) + 4g k T p Ju k ,v k /k 2 ,w k /k 4 
+ k 2( - n ~ 2p+ V [dfg k f p n>0fi (u k , v k /k 2 , w k /k A ) + 2g k T p n;0 (u k , v k /k 2 , w k /k A 

All numerical evaluations require a particular choice of the cutoff-function R(°'(z). 
the purpose of this paper, we will restrict ourselves to the use of the optimized cutoff 



(C6) 



or 



37 1 



where R, 



(o) 

opt 



(1 — z)9(l — z). The main virtue of this choice of cutoff is, that all the 



integrals appearing in the (generalized) threshold functions can be carried out analytically. 
In particular, 

1 1 =„, , 1 1 



r(n + 1) (1 + w)p ' nX ' T(n + 2)(l + w)P 
Similarly, the generalized threshold functions (IC4I) become 

1 1 
1 in + 1) (w + w + iy ) 



^ = ^ln(fl fc )=2. (C7) 



;-l) n + + + !,/ + !) 



(C8) 



r(n) 



(it + u + w) p 
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Here, /?(— l,n, I) denotes the incomplete beta function. For fixed values n, I, these become 
constants. This property leads to considerable simplifications in the analysis of the corre- 
sponding /3-functions. 



2. Evaluation of the traces 

The evaluation of the operator traces proceeds by expanding the arguments in a Taylor 
series in R around R = 0, keeping terms up to R 2 only. The operator traces appearing as 
"expansion coefficients" can then be evaluated with the heat-kernel techniques introduced in 
the last subsection, c.f. eq. (jCip . In particular, the functionals Q n [W] arising in these cases 
are of the form ( 1C5I) or ( 1C6I) . so that expressing them in terms of the generalized threshold 
functions is rather straightforward. We will now give the results for this evaluation for the 
various traces appearing in the main part of the paper, projecting the resulting RG flow 
onto the subspaces spanned by our truncations. 

The functions Q ra [W] featuring in the universal traces Sq and «Sit and the matter trace 
( TT8|) are a special case of (105j) with = and gj~ a ^-independent constant, which we can 
set to one. Following the strategy outlined above, the evaluation of the matter trace i fTBl 
results in 



^matter — t A \ 2 / ^ X \[9 



k A <$>\ + -k 2 $\R + ^(— R uvptT W a ' pa + — R 2 ) 
2 6 1 2 v 180 M p 80 ' 



(C9) 



(4vr 

while the expansion of the universal traces (|T9|) results in 
So = -^J d 4 *^ [k^ 2 + §(*} + 2& 2 )k 2 R + + + \<S>\)R 2 + ^R pvpa R^] , 

Sit = - J d A x^ [3k% l 2 + <$>\k 2 R + J|^ 2 - %R pvpa R^ p °] . 

(CIO) 

Here, all the $^ are evaluated at zero argument, and for the optimized cutoff they are 
trivially obtained from flOT[) . 

The non-universal traces entering into the computation of the /3-functions in the Einstein- 
Hilbert truncation (121!) can be evaluated along the same lines, utilizing the general relation 
(IC5j) . The results then read 

S hh = / d*x^g\k* fe(-4A/3) + ^ln(« 1 )$i(-4A/3)) 

W J / (en) 

+ |& 2 i2 ($}(-4A/3) + i9 t ln( Ml )<li(-4A/3) N " 
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and 



^ / d^hk 4 fe(-2A) + ldM Ul )$l(-2\) 

(4tt) 2 y l v / ^ (C12) 

- 5k 2 R (±$}(-2A) + |^1(-2A) + \d t \n{ Ul ) (|$|(-2A) + §$|(-2A;f 



respectively. 

The evaluation of the non-universal traces (|26|) appearing in the R 2 + C 2 truncation, 
on the other hand, is slightly more involved. Applying flC6p . the expansion of Shh on the 
truncation subspace takes the form 

Shh = 2{Ly / ^Va^^ + ^eR+^CsR^R^ + C.R 2 ^ . (C13) 

The dimensionless coefficients C, can be readily expressed in terms of the general- 
ized threshold functions ( 1C4I) with arguments T£ m = T£ jm (6<7 2 , #1, 2/3go), T p nml = 
^n,m/%2, 9i,2/3g ) and read 

d = 24<7 2 T 2;1 + 2 5l T^ + Qg 2 tl 0A + (2 9l + 9l )f * 0)0 , 

C 2 = 4<? 2 (l2y 2 T 2 2j2 - T 2j0 + + ^T 2 , x ) + i^T^o 

- g 2 (2f l A0 - ^Tjj^ - f } A1 ) + (2^i + <7i)(2<7 2 T 2 1>0 + |T^ 0j0 ) , 
(12^2 + ^i)^ + 6^2 + 2^i + ^i 

6^2 + Pi + 3^0 

C 4 = g 2 {96g 2 2 rl. 3 - S^T 2 , x - T 2 _ 2 - 5l T| 2 ) + f^T 2 , - T} >0 )} 
+ p 2 {24p 2 t| 2il - 2p 2 (2f 2 ilj0 - f 2 u ) - |f } A0 } 
+ (2^i + <h) (4<fef | 2)0 + §T 2 1(0 ) + ^C 3 . 

Here, we applied i?^ ^ (0) = 1 to simplify C 3 and expressed the resulting coefficients in terms 
of the dimensionless coupling constants (128]) . 

The projection of S 2 t proceeds in a similar fashion. In this case all threshold functions 
appear with arguments T* m = TP >m (2g 3 , -g u -g ), f p n m l = f p n m l (2g 3 , -g u -g ). Param- 
eterizing 

S 2T = j d^^Ct + C 2 k 2 R + ^-C 3 R aP , u R a ^ v + C 4 i? 2 | , (C15) 
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the coefficients Cj appearing in the trace expansion are 

C 1 = 4Qg 3 Tl A - 10^T^ + 10<j 3 f * 0>1 - 5(2 9l + gifil,o,o , 

C 2 =10^(4^ 3 T| 2 - - g l Tl 1 ) + 5^(4^^ + §T* - 9l T 2 2>0 ) - f g 3 T\ tl 

+ 5^ 3 (2^T| lil + g{T l fi>1 - |T} 01 ) - 5(2 gi + gi)(g^f 2 2>lfi - |f J 00 + \g{T 2 00 ) 

p? _ ( 4 ^3 - gi)(p + 2^3 - (2gi + gi) 

°3 — o ' 

2# 3 -gi-go 

C A =5g b9l {8g 3 Tl 2 - TC* - 2g 1 T% 1 + fT^} + |V {T} )0 - 4<? 3 T 2 )2 } - f 
+ 5g 2 {8g 3 T 3 2j3 - 2T\ X - 2 9l Tl 2 } + 5gf {2g 3 T 3 2jl + |T? >0 - 1^T» } 
+ 5 (2 5l + [g, (if ? >1)0 - yi f 1 )1(0 - (fcf | 2)0 ) + |^if 2 1AQ - i^T| 0)0 } 
+ 5<? 3 {^(2^T| 2jl + 2 5l f l hl - |f - | 5 if 2 , ,i + |^i 2 T| 0)1 } 

- 5^b jfl'bf 2,1,0 + Ifl'l ^2,0,0 ~~ I^IAOJ _ if^ • 

(C16) 

Note that the generalized threshold functions entering into C{ and C{ depend on different 
arguments. 

Substituting these expressions into the generic form of the flow equation (JTTJ) and com- 
paring the coefficients on the left and the right-hand-side then gives rise to the /3-functions 
(j22p for the Einstein-Hilbert case and fT29l) for the C 2 + /2 2 -truncation, respectively. 
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